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Abstract 

We present a general approach to study the flooding time (a measure of how fast infor- 
mation spreads) in dynamic graphs (graphs whose topology changes with time according 
to a random process). We consider arbitrary ergodic Markovian dynamic graph process, 
that is, processes in which the topology of the graph at time t depends only on its topol- 
ogy at time t — 1 and which have a unique stationary distribution. The most well studied 
models of dynamic graphs are all Markovian and ergodic. 

Under general conditions, we bound the flooding time in terms of the mixing time of the 
dynamic graph process. We recover, as special cases of our result, bounds on the flooding 
time for the random trip model and the random path models; previous analysis techniques 
provided bounds only in restricted settings for such models. Our result also provides the 
first bound for the random waypoint model (which is tight for the most realistic ranges of 
network parameters) whose analysis had been an important open question. 
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1 Introduction 



A dynamic graph is a probabilistic process that describes a graph whose topology changes with 
time. Dynamic graphs are appropriate models of wireless networks, peer-to-peer networks, 
social networks, and so on. There are several interesting problems on dynamic graph processes, 
for example load balancing, studied in [161 EH]- Here, we are interested in the speed of 
information spreading, a question that can model the spread of disease, the broadcast of files 
on peer-to-peer networks, of memes in social networks, etc. 

The simplest model of information spreading is the process of flooding [H [lOl [131 13 [13 
[TSl [22l [26l I27j . which begins with one node in the network being given a certain piece of 
information, and then the application of a protocol in which, at each time step, every node 
that has the information spreads it to its neighbors. Recall that the neighborhood of a node 
changes with time, and so even though the flooding algorithm is deterministic, the process of 
information spread is probabilistic. 

The setting of wireless networks has motivated the study of this problem in geometric 
models of dynamic graphs in which, at every time step, every node is mapped to a point 
in a metric space, and two vertices are connected if their distance is smaller than a given 
communication radius. The underlying metric space is usually a bounded portion of the 
plane, for example a square or unit disk, and the dynamics come from independent random 
walks performed by the individual nodes via local moves [201 [13 IHl [121 [Ml [23 [23]. For 
example, a representative model of this type is the random walk model: n nodes are placed 
on a m X m grid; at each time step, every node v independently moves to a point in the grid 
randomly chosen among the points adjacent to the one that v occupied at the previous time 
step; at each time step, the edge {u, v) is present in the dynamic graph if u and v are located 
within distance r in the grid. In such a model, the flooding time will depend on the initial 
locations, on r, m and n. Usually, one is interested in a worst-case analysis with respect to 
the initial locations, and in a bound dependent only on r, m and n. Known analyses of such 
models rely rather strongly on the fact that, if we consider a fixed node v, the stationary 
distribution of its location (i.e. the positional distribution) on the grid under the random 
walk is essentially uniform. 

The random waypoint model [3 [3 [3 [21] is another classic model of networks of mobile 
agents, and it is probably the most well studied one. In this model, every node chooses a 
random destination point in the mobility space, then he travels over the shortest path till 
he reaches the destination, and so on. Some analytical properties of this model have been 
derived such as the mixing time, the stationary node distribution, etc. [3 [3 121] • However, 
bounding its flooding time (or any basic communication tasks such as routing, data collection, 
etc.) is still a fundamental open problem. The techniques adopted for the random walk model 
do not work in the random waypoint model, mainly because of the presence of long periods 
the node spends in deterministic trajectories and because of the strong difference between 
their respective stationary positional distributions. For instance, the stationary positional 
distribution of the random waypoint over a square is in fact far from the uniform one: it is 
highly biased towards the center of the square. 

The above two families of models define a probabilistic process over the nodes, which then 
implies which pairs of nodes are connected by an edge. There are also models of probabilistic 
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processes that are directly defined over the edges. A very general model is provided by 
Markovian Evolving Graphs (MEGs) introduced in [2] (see later for a formal definition). 
However, available techniques to analyze information spreading only concern very restricted 
subclasses of MEGs, such as that studied in [10], where the state (i.e. on/off) of the edges is 
ruled by independent copies of a simple two-state Markov chain. 



Our Work 

- General Dynamic Networks. We provide an upper bound to the flooding time of any sta- 
tionary Markovian Evolving Graph (MEG) [2]. If we call Gt = {V,Et) the random variable 
describing the dynamic graph at time t, the process is a MEG if the sequence of random 
variables ^([n], {Et}t^o) = Gq, . . . , Gt, ... is Markovian, that is, if the distribution Gt is com- 
pletely determined by the distribution Gt-i, via a transformation independent of t. The class 
MEG is extremely general and, in particular, it includes all above-mentioned network models 
such as random walk and random waypoint: indeed, it is easy to verify that any of such mod- 
els yields a sequence of random graphs G{[n],{Et}t^o) which is Markovian. We require the 
Markovian sequence to converge to a unique stationary distribution, for every initial choice 
Gq. This stationary distribution is a probability distribution over n- vertex graphs and it will 
be called the stationary graph Qn- 

Our upper bound to the flooding time is a function of: i) the edge-probability in Qn (that 
determines the expected density of the stationary graph) ii) the degree of independence among 
ed ges in Qn and, iii) the mixing time of the Markov chain ^ ([??,], {^/^Ij^g). More precisely, 
given edge e, node i, and node subset A, let p{e) be the probability that e exists in Qn and 
let e{i, A) be the binary random variable returning 1 iff an edge exists connecting i to some 
node in A in Gn- Let M be an upper bound on the mixing time of ^([n], {Et}t^o)- Our result 
states that if, for some arbitrary positive reals a and /?, (i) p{e) ^ a for any link e anc|l| 
(ii) P (e(i, A) • e{j,A)) ^ (3P (e{i, A)) ■ P (e{j,A)) for arbitrary nodes i,j and arbitrary node 
subset A (not containing i or j), then the flooding time is w.h.pH 

Our methods can be applied to non-Markovian processes as well, although the results are 
more complex to state (see Sections [21 [3|) . Bounding the mixing time of dynamic networks has 
been the subject of several studies in the last years [Il[24l[in]: our result allows to efficiently 
exploit any (previous or future) bound on the mixing time for the flooding time of such 
MEGs. In order to get an intuition of the real meaning of Conditions (i) and (ii) (i.e. how 
mild they can be), we observe that mild bounds on the density and independence parameters, 
say a = Q{l/n) and (3 = O(polylogn), do not imply any good node/edge expansion of the 
single snapshot graphs Gt of the process: In every Gt there could be a large subset of all nodes 
(say half of them) that are isolated. In such sparse and disconnected topologies, thanks to our 
bound, the flooding time can be just a poly-logarithmic factor away from the mixing time of 



^With an abuse of notation, event probabilities such as P (e(i, A) = 1) will be shortly denoted as P (e(i, A)). 
^We say that an event holds with high probability if it holds with probability at least 1 — 1/n. 
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the MEG. This crucial fact has strong consequences on concrete network scenarios which are 
discussed below. 

- Node-MEGs. Our general approach finds a natural application in a subclass of MEGs that 
we call node-Markovian Evolving Graphs, denoted as node-MEGs, where every node changes 
its state independently according to a Markov Chain Ai. The state of a node can implement 
several dynamic features of the node (such as geometric position and destination, trajectory 
phase, social role, etc). Then the existence of a link between two nodes (only) depends on 
the current states of the two nodes according to a fixed deterministic function. Observe that 
node-MEGs are a relevant class of MEGs that includes every mobility model where nodes acts 
independently over any discrete space (such as an arbitrary graph or any geometric region 
discretized, for instance, by using a grid of suitable resolutioro) . Random walk, random 
waypoint, and random trip models [71 [231 [H] have a natural realization as node-MEGs: for 
instance, a realization of the random waypoint as a node-MEGs is described in Subsection 14. 11 
Notice that in node-MEGs, nodes are indistinguishable so Condition (i) is easy-to-check: if it 
is satisfied (in the stationary graph) for a specific edge then it is satisfied for all edges. We then 
prove that if incident edges are almost pairwise independent in the stationary graph yielded 
by a node-MEG then Condition (ii) is satisfied for some constant f3: so, checking Condition 
(ii) can here be reduced to check pairwise independence of incident edges in a generic node. 

- Geometric Mobility Models. As for node-MEGs defined over geometric spaces (such as the 
random waypoint), we state an easy-to-check condition implying Condition (ii): we transform 
the pairwise-independence condition on incident edges into some mild uniformity conditions 
on the single-node positional stationary distribution (see Corollary |4|) . We in fact show the 
former properties are satisfied by a wide class of random mobility models such as the random 
waypoint over a square: we thus get the first known upper bound for its flooding time. 
Furthermore, when the stationary graph is sparse, the obtained bound is almost tight, i.e., 
0((L/f )polylogn) where L is the diameter of the square and v is the node speed. In particular, 
our bound is almost tight whenever the transmission radius and the node speed are absolute 
constants: this is surely the model setting that best fits opportunistic delay-tolerant Mobile 
Ad-hoc Networks [Ml [fsi ET] . 

- Graph Mobility Models. Our upper bound holds even when the mobility space is an arbitrary 
graph (the vertices of such graph are called points): in this case nodes choose randomly their 
trajectories from some fixed families of simple paths of the mobility graph and the mixing time 
is proportional to its hop-diameter. This mobility model is called random paths (on graphs). 
The parameter /3 in Eq. [T]is here determined by the point congestion yielded by the feasible 
paths of the mobility graph: informally speaking, /3 is small whenever the random paths do 
not yield a high point congestion (again, we get a somewhat mild uniformity condition on the 
(single-node) positional stationary distribution). If this is the case, the obtained bound on 
the flooding time becomes almost tight. The random walk model on graphs can be seen as 
a special case of the random paths on graphs: then the obtained bound on the flooding time 
improves over the previous one [15J for the class of graphs where the mixing time of a random 
walk is shorter than the meeting time of two random walks. 

^The level of resolution does not afltect the obtained bound on the flooding time, provided the resolution is 
high enough. 
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By concluding about mobile networks, we want to emphasize the impact of our general 
method over classic models: the general conditions (i) and (ii) of our bound are transformed 
into mild uniformity conditions over the positional stationary distribution which can be verified 
by standard techniques [U [23] . 

- Link-based Dynamic Networks. In appendix we show that our method can also be applied 
to a subclass of MEGs in which each edge evolves independently according to an arbitrary 
(hidden) Markov chain. Previously \10\ [5]. such link-based dynamic models had been studied 
only in the case in which the edge Markov chain is very simple. 

Previous Works. As mentioned before, information spreading in dynamic networks has been 
extensively studied in the literature under a variety of scenarios and objectives (for a recent 
good survey see [22] )• For brevity's sake, we restrict our attention to the results more directly 
related to our work. Previous models and results can be roughly classified in two main classes: 
link-based dynamic graphs and mobility models. 

As for the first class, in [9], radio broadcasting is analyzed on a dynamic graph managed by 
a worst-case dynamic adversary and on a sequence of independent Erdos-Reny graphs. In [10] . 
an upper bound on the flooding time for the restricted model edge-MEG has been derived. 
The flooding time of another simple version of edge- MEGs has been studied in [5|. The general 
MEG model has been introduced in p] where some results are obtained for the cover time and 
hitting time of random walks. Flooding time in stationary MEGs is studied in [11]; unlike our 
method, the method in [11] only holds for stationary graphs which are connected and good 
expanders. A worst-case model of dynamic graphs has been introduced in [21]. The analysis of 
some communication tasks is presented under the strong stability condition called T-interval 
connectivity (for T ^ 1) which stipulates that for every T consecutive steps a stable connected 
spanning subgraph must exist. 

As for mobility models, almost tight bounds on the flooding time for the random walk 
model have been obtained in [201 El 121 [Ml [271 US]- As above discussed, their techniques 
strongly rely on specific properties of the adopted version of the random walk. The case of 
general mobility graphs has been considered in [15]: the obtained results are discussed and 
compared to our results in Subsection 14. 1[ An upper bound on a variant of the random 
waypoint model has been derived in p^. In this version, nodes follows Manhattan paths. 
Similarly to the works on the random walk models, the ad-hoc analysis in [13j strongly relies 
on the particular node trajectories and on the specific positional distribution yielded by this 
model. So, its contribution strongly departs from our general approach that obtains bounds 
for any version of the random waypoint model. 

Paper Organization. Section [2] formalizes the general model of dynamic graph and the 
flooding process. Section [3] provides the main theorem for the flooding time in the general 
model. The node- MEG model is described in Section [Jj where an upper bound on the flooding 
time is given for this specific model. The representation of the random trip and the random 
paths models as specific instances of the node-MEG model is given in Subsection 14.11 here, 
the flooding-time bound for node-MEGs is transformed into two useful bounds on the flooding 
time: the first one for the random trip and the second one for the random paths. Finally, 
conclusions with open questions are discussed in Section[5l Due to lack of space, the application 
of the main theorem to general edge-MEGs and all the proofs are given in the appendix. 
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2 Preliminaries 



In this section, we introduce the general model of dynamic graphs. For any positive n, [n] will 
denote the set {1, 2, . . . , n}. A dynamic graph G{[n], {Ef}f^Q), with node set [n], is a stochastic 
process represented by a sequence of random variables Eq, Ei, . . . , Et, . . . such that, for every 
t, Et is the set of edges of the dynamic graph at time t. The speed of information spreading 
can be studied in terms of the flooding time. Given a dynamic graph Q{[n],{Et}t^o) and a 
node s € [n], the flooding process with source s is defined as follows. At time t = 0, ,s is 
the only informed node; then a node v gets informed at time t + 1 iff an edge e £ Ef exists 
connecting some informed node to v. Flooding over a dynamic graph is represented by the 
stochastic process {It}t^o where 

lo = {s} and yt^llt = It-i U {j e [n] \ 3i G h-i : {^, j} G Et} 

The random variable It is the set of informed nodes at time t. Clearly, it holds that 

/o ^ /i c /2 c . . . It c . . . 

The flooding time with source s is the random variable F{Q,s) = min^j/t = [n]} and the 
flooding time is the random variable F{Q) = max^ F{Q, s). 

Given a dynamic graph G{[n\,{Et}t^o), we define the following random variables. For every 
time t, for every pair of nodes i,j G [n] and for every subset of nodes AC [n], let 

^t^( 1 if{i,i}eEt e* =1 ^ ii3jeA:{i,j}eEt 

''^ \ otherwise \ otherwise 

Moreover, for any binary random variable X, the notation P (-'^ | Et^r) ^ (or ^) a stands 
for 

V sequence of edge subsets Aq, . . . , At, with P I /\{Et = At) ] > 0), it holds 

PIX = 1 /\{Et=At)] <(or^)a 
t=o / 



3 Flooding in Dynamic Graphs 

Our goal is to evaluate the flooding time of a dynamic graph as a function of some properties 
of its edges. These properties are not required to show up at every snapshot, rather it suffices 
that they hold at the beginning of every time "epoch", where an epoch is a sequence of 
consecutive time steps of suitable length. When the dynamic graph is a Markovian process 
admitting a stationary graph, the properties above refer to the expansion properties of the 
stationary graph and the epoch length refers to the mixing time. However, aiming at the 
maximal generality, we introduce such concepts for general (non-Markovian) process. 

Let M be a positive integer and let a and /3 be two positive reals. A dynamic graph 
Q{[n],{Et}t^o) is (M , a, P)- stationary if Vr ^ 1, Vi,j G [n] with i / j, \/A C [n] — the 
following two conditions hold 
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1. P (e[f I ^;t^(r-i)M) ^ a (Density Condition) 

2. P (e[_f • ej*^ I ^(^(,_i)A/) ^ I3P (e[l | ^t^(,-i)A/) P (ej^^^ | ^t^(.-i)M) (/^-Independence 
Condition) 

The full proof of the following theorem is given in the Appendix. 

Theorem 1 (Flooding Time) IfQ{[n],{Et}t^o) is {M , a, f3)-stationary then with high prob- 
ability the flooding time in Q is 

O (^-^ + log2 

4 Node Markovian Evolving Graphs 

We introduce the general class of Node Markovian Evolving Graphs (in short, node-MEGs) 
where the behavior of the nodes is ruled by independent Markov chains. To every node is 
associated a Markov chain whose states contain enough information to determine whether two 
nodes are connected or not. Of course, this is an approach based upon hidden Markov chains 
to model dynamic graphs. 

Let M. = {S, P) be the Markov chain associated to every node, where S is the set of 
states and P : S x S ^ M are the transition probabilities. The connections are determined 
by a symmetric map C : 5 x 5 — )• {0, 1}: any two nodes i, j are connected at a given 
time t if C{u,v) = 1, where u, v are the states of i and j at time t. The symmetric map 
C(-, •) is also called the connection graph of Ai. A node-MEG is denoted by NM{n,Ai,C). 
Notice that the Markov chain M may depend on the number of nodes. The initial state of 
each node i is random with a probability distribution ^ over the set of states S. We denote 
by L the global initial probability distribution determined by the product of the probability 
distributions Li. The state of a node i at time t is a random variable s* fully determined by the 
initial distribution and the Markov chain M. A node-MEG NM{n,M,C) together with 
an initial probability distribution l determines a dynamic graph ^([ra], {Et}t^o) where, for any 
t 0, Et = {{i,j} I C(s*,s*) = 1} It is easy to verify that node-MEGs enjoy the following 
property. 

Fact 2 Consider any node-MEG NM = NA'I{n,A4,C) in its stationary state, then the prob- 
ability -Pnm that any fixed pair of nodes are connected and the probability PnM2 that two fixed 
nodes are both connected to another fixed one do not depend on the choice of the fixed nodes: 
they are functions only of the stationary distribution of M and of the symmetric map C(-, •). 

In Subsection 14.11 we will show that a wide class of mobility models turns out to be a special 
instance of Node-MEGs. 

Flooding in Node-MEGs. We now derive some simple properties ensuring that a node- 
MEG NM{n,A4,C) is a (M, a, /?)-stationary dynamic graph. Since in a node-MEG, edges 
are not independent, the crucial condition is the /3-independence. The models at hand are 
Markovian so the idea is to consider the model during its stationary state, that is, the time M 
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between two consecutive epochs is proportional to the mixing time of the Markov chain M. 
Moreover, the /3-independence involves sets of incident edges of arbitrary size (i.e. the size of 
subset A): instead, thanks to the independence of node evolutions and Fact [21 in node-MEGs 
this condition can be relaxed to a parameterized pairwise independence among incident edges. 

Theorem 3 Let NM = NM{n,M, C) be a node-MEG such that Pnm ^ l/n^^^) and Pnm2 ^ 
^(-Pnm)^; for some rj ^ 1. Then, with high probability, the flooding time is 



where Tmix is the mixing time of the Markov chain Ai . 
4.1 Flooding in Classic Mobility Models 

Geometric Mobility Models. Several mobility models can be represented as special cases 
of node-MEGs. Many of these are continuous-space models [71 [23] in which nodes move 
over a subset of W^. Since node-MEGs are discrete, we approximate continuous space by 
discretization. In the simplest and most common case, nodes move over a square of of side 
length L. The square can be discretized by taking a square grid Q formed by m x m points 
regularly spaced in the square region, where m can be arbitrarily chosen. 

In the standard random waypoint [7], n nodes independently move over the square: every 
node randomly chooses a speed in [vmin-,Vmax\ where Vmax = ^{vmin) and a destination point 
('waypoint') in the square and moves with the chosen speed on a straight path to this point. 
Then, it repeats the same process again and again. The destination points are uniformly 
distributed on the square. At any time two nodes are connected if they are at distance not 
larger than the transmission radius r. 

The formal discretization of the random waypoint as a node-MEG can be done by simple and 
standard arguments [H], so we here provide an informal description only. The generic state 
of the Markov chain M must encode the destination point, the current point in the straight 
point-path the node lies, and the node speed (the latter can be defined as the number of 
points per time step). Then the transition matrix can be easily defined: when a node is in 
some internal point of a path the choice of his next state is deterministic while when he arrives 
at the end of a path, his next state is randomly chosen by selecting the next destination point 
(and thus the next path to be followed) and the speed. As usual, the connection map C is 
defined as follows: there is an edge between nodes u and v at time t iff their relative distance 
at time t is not larger than r. The positional probability distribution in the stationary phase 
is defined as the probability that a node is in point x (for any choice of x over the square) 
when the state of the node is random with the stationary distribution of M. The density 
function of this distribution (in short, positional function) yielded by the random waypoint 
over the square will be denoted by Fyjp{-). The random waypoint belongs to a general class 
of geometric mobility models called the random trip model [24] where the mobility space IZ 
can be any bounded connected region of and the feasible node-trajectories can be any 
family of continuous curves. Any random trip model can be discretized with any level of 
"approximation" (in terms of grid resolution and time unit) by following the same procedure 
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described above for the random waypoint. For this geometric class of node-MEGs, Theorem 
[3] can be rewritten in a very useful way. 

For any r ^ 0, D{u, r) denotes the set of all the points that are at Euclidean distance at most 
r from u. For any connected region B C M*^, define Br = {u ^ B \ D{u, r) C B^ and vol(;B) 
as the volume of region B. The proof of the next bound is given in the Appendix. 

Corollary 4 Let NM = NM{n,Ai,C) be a node-MEG yielded by a suitable discretization of 
a random trip model T over a bounded connected region TZ QW^ with positional function Ff. 
If Pnm ^ l/n*^^"*^) and for some 6^1 and X > it holds that 

(a) Vn e n, Fr{u) ^ ^ 

(b) 3Bcn such that vol{Br) ^ Avol(7^) and Vn G B, Fr{u) ^ j^;^ 

then, with high probability, the flooding time is O ^mix + p') log^ where Tj^ix 

is the mixing time of the Markov chain M . 

The useful novelty of the above corollary lies in the following fact: the pairwise-independence 
condition in Theorem [3] is transformed into two mild "uniformity" conditions on the positional 
function yielded by the mobility model. The latters only refer to the stationary positional 
distribution of the single node and it is often much easier to verify with respect to the pairwise 
condition. Indeed, a general method (the Palm Calculus) to derive explicit formulas of such 
function for random trip models has been introduced in [24]. As for the random waypoint on 
the square, the explicit positional function -F^p(-) has been derived in [25j and it is easy to 
verify that the two conditions of the above corollary are satisfied for some absolute constants 
5 and A. Furthermore, the mixing time of the random waypoint over a square of side length L 
is Q{L/vmax) (remind we are assuming Vmax = 0{vmin)) [IlIlH]. We thus obtain the following 

bound on the flooding time O f — ^ + 1^ log^n | Let us consider the case L ~ -v/n, 

r = 0(1), and r = 0{v 

max)'} notice that this standard setting yields a stationary mobile 
network which is w.h.p. sparse and highly disconnected. Then the bound on the flooding time 
becomes O ( ^ log^ n ] which almost matches the trivial lower bound 57 ( ^ ) . 

Graph Mobility Models. A natural generalization of random walks over a graph can be 
defined by considering random paths over a graph. This clearly includes the random waypoint 
over a graph. At every time step, a node moves along a path instead of on a single edge. More 
precisely, the model is specified by a graph Hiy, A) and a family V of feasible paths in H 
satisfying the property: for every path h £ V, there is a path h' G V such that h' starts where 
h ends. For any u € V, let V{u) be the set of paths in V that starts at point u. The mobility 
model is as follows, a node at point u G V chooses uniformly at random a path in V{u), then 
it travels along the path (an edge at the time), when it reaches the end point v, it chooses 
uniformly at random a path in V{v) and travels along that path, and so on. We assume that 
two nodes are connected, at any given time t, if they are in the same point at time t. For 
any path h, let i{h) denote the number of points of h. The representation of a random path 
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model RP = {H,V) as a node-MEG is straightforward. The Markov chain A^rp = {S,P) is 
such that S = {{h,hi) \ h£V,2^i^ l{h) and hi is the ith point of /i}; the transition 
probabihties are as follows 

V/i e 7^ Vi : 2 s$ i < (ih) P{{h, hi), {h, hi+i)) = 1 

1 



V/i, h' eV: ht(h) = K P{{h, heih)), {h\ h'^)) = |^(/^^^^^)| 
all other transition probabilities are equal to zero; 

and the connection map is such that Crp((/i, /ij), (/i', h'-)) = 1 iff /i, = h'-. Observe that if V is 
the set of edges of H then the mobility model is equivalent to the random walk over H (with 
p = l). 

We say that a path h £ V passes through a point u \i hi = u for some 2 ^ i ^ l!.{h). For 
any point u &V, let i^-p{u) be the number of paths in V that passes through point u. Notice 
that if V is the set of edges of H, then #-p(u) = degjj{u). A random-path model RP = {H, V) 
is said to be simple if every path in V does not pass through the same point more than once, 
but the start and end points that may be equals. Moreover, RP is reversible if, for every path 
h £ V, the reverse path of h also belongs to V. We say that a random path model {H, V) is 
5-regular if 

V-u G F #v[u) ^ 5 1^-^ 

Roughly speaking 5-regularity ensures that no point is a much busier crossroad than the aver- 
age. Theorem [3] then implies the following useful result (the proof is given in the Appendix). 

Corollary 5 Let NM = NM{n, Ai, C) be a node-MEG yielded by a random path model RP = 
(H,V) that is simple, reversible, 6-reqular, and \ V\ ^ Then, w.h.p the flooding time is 

o(.,.(e.,,.„..)..„.„„..... 

If for every pair of points there is only one feasible simple path, then the mixing time of the 
relative Markov chain is O(-D), where D is the diameter of H. Moreover, if the model is 
(5-regular for some 6 = polylog(n) and \V\ = 0(npolylog(n)), then the above corollary implies 
that the flooding time is O {D polylog(n)). This is within a poly-logarithmic factor from the 
optimum, since a trivial lower bound is Q(D). A basic instance of this case is when H is a 
grid and the feasible paths are the shortest ones. 

As mentioned in the Introduction, almost tight bounds for the flooding time on the random 
walk model over grids have been recently obtained in [HI [121 [27] . In ^ general graph H(y, A), 
every node randomly chooses his next position among all points in V that are within p hops 
from his current position. The transmission radius r determines the maximal distance (again in 
terms of number of hops in H(y, A)) within which a message can be successfully transmitted. 
The most studied setting is p = 1 and r = 0: a node makes at most one hop per time step and 
it can infect only nodes that lies in the same point. This natural setting in general graphs has 
been studied in [15]: the flooding time is proved to be 0(T* log n) where T* is the meeting 
time between two independent random walks on H. 
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In what follows we apply our analysis for the random path model to the special case of 
random walks. The (5-regularity condition over paths is transformed into a simple condition 
on the degree of the points. Given any 5 ^ 1, a graph H(y,A) is said to be (5-regular if 
(max{deg(t>) \v G y})/(min{deg(t;) \v £ V}) ^ 6. Then, we can derive a simple adaptation 
of Corollary [5l 

Corollary 6 Let NM = NM{n,Ai, C) be a node-MEG yielded by the random walk over any 5- 
where Tj^m is the mixing time of a random walk over H . 

The above bound improves the result in [15] over a relevant and wide class of mobility graphs. 
Indeed, given a symmetric graph, the meeting time of two random walks is asymptotically 
equivalent to the first hitting time and can be much larger than the mixing time of a single 
random walk [1]. A natural example is that of /c-augmented grids: take a grid of s points and 
add an edge between any pair of points whose hop-distance is not larger than k. While the 
meeting time is not smaller than that of a standard grid (slogs) [U [27], the mixing time 
decreases in k. For instance, if s ~ repolylog(n) then the bound in [15] becomes 0(npolylog(n)) 
while our bound is 0((npolylog(n))/A;^). 



5 Conclusions 

We believe that significant improvements are possible to the bound in Theorem [Tj along some 
interesting directions. We suspect that, under mild assumptions on the dynamic graph process, 
the factor + can be improved. A more challenging task is to avoid the dependency 
of the bound on the mixing time of the graph process. The density and /3-independence 
conditions can be met even at a state in which the graph process is far from the stationary 
distribution, and so a more refined analysis might be able to bound the flooding time without 
having to "wait" for the process to achieve stationarity. 

Our method may prove useful in the analysis of more refined communication protocols than 
flooding. A simple example is a randomized protocol in which, at every step, a node that 
possesses the information transmits it to a randomly chosen subset of neighbors. The analysis 
of such a process can be reduced to the analysis of flooding in a "virtual" dynamic graph in 
which a subset of the edges are removed. More general communication protocols might also 
be reduced to flooding by folding the actions of the protocol into the dynamic graph process. 
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A Generalized Edge-MEGs 



The link-based dynamic model Edge-MEG has been introduced in |10j and successively studied 
in [HIIIlIITj. In this restricted instance of MEGs, edges evolves independently. Every edge of 
the n-node graph can be in two states only: on or off. At any time step, every edge changes 
its state according to a two-state Markovian process with probabilities p (edge birth-rate) and 
q (edge death-rate): if at time t an edge e exists then it will die with probability q while if e 
does not exists then it will come up with proability p. In [T0|, the authors prove an almost 
tight upper bound on the flooding time. 



A more refined model with four states has been recently introduced and studied in [5]. 
Edge-MEGs do not explicitly model node's mobility, rather they are more suitable to model 
the link evolution in peer-to-peer networks or faulty networks. 

Our main contribution here lies in the fact that Theorem [1] can be applied to the much 
more general version of edge- MEGs where an arbitrary (hidden) Markov chain Ai = {S,P) 
rules the behavior of every edge and by an arbitrary map x : 5 — t- {0, 1} that determines, in 
function of the state, whether the edge exists or not. This wide generalization of edge-MEG 
will be denoted as EM{n, Ai,x)- 

The initial state of each edge {i,j} is random with a probability distribution i^ijy over the 
set of states S. We denote by l the global initial probability distribution determined by the 
product of the probability distributions '■{ij}. The state of an edge {i,j} at any time t is a 
random variable sj- completely determined by the initial distribution and the Markov 
chain Ai. 

Any model EM{n,Ai,x) together with an initial probability distribution l determines a dy- 
namic graph ^([ra], {Et}t^o) where, for any t ^ 0, 



A crucial property of such generalized edge-MEGs is that edges are independent random 
variables, so it always holds that the /3-indpendence is satisfied with /3 = 1. Then, when the 
Markov chain Ai admits a unique stationary distribution tt. Theorem [T] implies that hooding 
time is 



where Tmix is the mixing time of the Markov chain Ai and a is the probability an edge exists 
in the stationary regime (i.e. according to vr). For instance, in the basic edge-MEG model 
with parameters p and q the mixing time is 0(l/(p + q)) and a = p/{p + q) [TO]. We thus get 
an upper bound 



By comparing our bound to the (almost-tight) one in Eq. O we get that the former is almost 
tight whenever q ^ np. 



0(log n/ log(l + np)) 



(2) 



Et = {{i,3]\x{sl,,}) = l] 
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B Some useful Inequalities 



The Paley-Zigmund inequality. If X ^ is a random variable with finite variance, and if 
< 6* < 1, then 

Pix;,eE[x]) ^ (3) 

Lemma 7 (Lemma 3.1 in [3j) Let Xi, . . . , X„ be a sequence of random variables with val- 
ues in an arbitrary domain, and let Yi, . . . ,Yn be a sequence of binary random variables, with 
the property that Yi = Yi{Xi, . . . , X„). // 

P{Y, = l\ 

then 

^(IZ^*^^) ^ F{B{n,p) i^k) 
where B{n,p) is binomially distributed random variable with parameters n and p. 

Lemma 8 (Chernoff Bound) Let be independent binary random variables and 

let X = Yll=i -^i ^'^d /U = E [X] . Then, for any 5 > 0, 

P (X < (1 - < exp 




C Proof of Theorem [T] 



Expansion Properties. In what fohows, we derive some expansion properties of (M,a,/3)-stationary 

dynamic graphs. Such properties will be then exploited in the analysis of flooding. 

The times tM will be called epochs and they will be abbreviated by r (i.e. r will stand for 

rM, with respect to a fixed (M, a, /3)-stationary dynamic graph ^([n], {E'jjt^o))- Thus, we 

write Er, ejj and to denote EtM, and e[^, respectively. These abbreviations will 

be also used for other random variables. 

Let degj^^ be the random variable counting the number of nodes in A connected to i at epoch 
r, i.e., deg[^ = \{j G A \ {i,j} £ ^r}|- Observe that, thanks to Condition (1), the expected 
value of degj^^ is at least \A\a; the following lemma provides a concentration result as function 
of the "independence" parameter /3. 

Lemma 9 LfQ{[n],{Et}t^o) is (Ad , a, f3)- stationary then, Vr ^ 1, Vi G [n], andMA £ [n] — {i}, 

P(degl,^^^ 



E<T--] 



\A\a 
2 + 2|A|a/3 
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Proof. Firstly we bound the expected square of the degree. For the sake of brevity, we omit 
the conditioning by E<^t-~i- It holds that 

v2l 



E 



(deg[,^)^ 



E 



j,keA 



-'1,3' 



= E[el,,.-ey +E[deg[,^] (since (e[,^.)' = < 

< E [el^] E [elk] + E [deg[,^] (by Condition (2)) 

^ /3E[deg[,^]' + E[deg[,^] 
Moreover, from Condition (1), it derives that 

E [deg[^^] ^ \A\a 

Thus, from the Paley-Zigmund inequality (with 6 = 1/2) and the above bounds, we obtain 

1 E[deg[,^]' 



P I deg[ 4 > 



\A\a 



2E[(deg[,^) 



^ -- 



E [deg- 



i,A\ 



2/5E [deg[,^]' + E[deg[,^] 



E [deg; 



Ia\ 



2 + 2/3E [deg[,^] 
E [deg[,^] 



2 + 2/3E [deg[,^] 
Since function 2+2px decreasing, we have that 



P I deg[ , ^ 



\A\a 



E [deg[,. 



\A\a 



2 ; ^ 2 + 2/3E[deg[,^] ^ 2 + 2/3|^| 



a 



n 

The next lemma extends Lemma [9] to the expansion of an arbitrary node subset. For every 
epoch T and for every A, i? C [n], define the random variable counting the expansion from A 
to B 

deg\B = \{j(^B\3ieA:{i,j]eEr]\ 

Lemma 10 If Q{[n],{Et}t^Q) is {M, a, /3)- stationary then, Vr ^ 1, \/A G [n], and S 
[n]-A, 



P deg^ o ^ 



\A\\B\a 
4 + A\A\af3 



E. 



\A\\B\a 



A + Q\A\\B\al5 
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Proof. The argument is very similar to that of the proof of Lemma [9l Firstly we bound the 
expected square of the expansion. For the sake of brevity, we omit the conditioning by E<^r-i- 



E [(deg^,^)' 



E 




j,keB 



^ E [degle] + /3 E E [4,^] ^ Ka] (by Condition (2)) 

j,keB 

= E[deg^,B] +/3E[degXB]' 



Moreover, from Lemma [9] it derives that 

E [deg^,^] = E [eU = P (degj,^ > O) ^ J] 



jeB jeB 
From the Paley-Zigmund inequality (with 



jeB 



\A\\B\a 
2 + 2\A\a/3 ~ 2 + 2|A|a/3 



\A\a 



l^ll^la 
'^^■^ ^ 4 + 4|A|a/3 



P deg^. R > 



1/2) and the above bounds, we have that 

1 E[degX^]^ 
2E[(deg^,^)2] 



E [deg; 



A,Bi 



2E [deg^^s] + 2/3E [deg^,^] 

E [deg^,^] 
2 + 2/3E [deg^^s] 

\A\\B\o^ 
2+2\A\oiP 

o l^l|g~ 

4 + 6|A||S|a;S 



n 

When the dynamic graph is sparse, the expansion rate obtained by considering a single snap- 
shot of the process (i.e. the expansion of a node subset at time r) does not suffice to get 
a good number of new informed nodes. In this case, a "dynamic" version of the expansion 
properties is required. For every epoch r, for every T ^ 1, and for every AC [n], define 

spread^^ = |{i G [n] - A | 3r'3i G A : r < r r + T A {i, j} e E^,}\ 

That is, spread^"^ is the number of nodes outside A that get connected to nodes in A during 
the epochs in the interval (r, r + T] . 
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Lemma 11 If Q{[n],{Et}t^o) is (M , a, (3) -stationary then, Vr ^ 1, Q [n] with \ A\ ^ n/i, 
and Vt ^ 0, 

P ^spread^'^ < \A\ E^r^ < exp(-t) 

where 

/3 , |A|/32^ 



T = 256 



1 



+ — + 



\A\na 



3/3 U 



Proof. For brevity's sake, we omit the conditioning by -E^r- Let St be the set of nodes outside 
A that get connected to nodes in A during the epochs in (r, t + t\. Formally, 

So = and St = U {j G [n] - A | 3i G A : {i, j} G 

Clearly, |S't| = spread^"^. Let 



7 



8 + 8|yl|a/3 



Define 



Yt 



1 if ^ 1^1 or ^ + [71 

otherwise 



Observe that = ft{Er+i, ■ ■ ■ , E^-^t) for a suitable function ft. From the inequality P {AV B \ H) ^ 
P ( S I H it derives that 

p(yt = 1 1 Er,t-i) = P(|5t-i| ^ 1^1 V \St\ > \St^i\ + h] I ^.,t-i) 

^ P(|5t| > \St-i\ + [71 I Er,t-i A < |A|) (4) 

where E^^t-i stands for E-^+i, . . . , Er+t-i- Assume that \St-i\ < \ A\ and let W = [n]—A—St-i. 
Since |^| ^ n/4, it holds that \W\ ^ ^ and 



\St\>\St-i\ + h] ^ deg^ti^^7 



From Lemma \T0\ we have that 



P(deg^+^^7 



E^^+t-i A < |yl| ^ P deg^+^ ^ 



|^||VF|a 
4 + 4|^|a/3 



(5) 



E^r+t-i/\\St-i\ < \A\ 



|^||l^|a 



4 + 6|^||VF|a^ 



^ + 6|^|na/3 
Thus, from Ineq.s|ll[5l and [6] we get 

P(Ft = 1 I Er,t-x) ^ P 
We can now apply Lemma [7] to the r.v. E'r+ii • • • > -E'r+T and r.v. Yi, . . . , 



(6) 
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Since spread^ < 1^1 =^ X^t=i < -rf--, from the above inequality we obtain 

P (spread^^ < ^ P {b{T,p) < 

By applying Chernoff's Bound (Lemma [8]) , after some calculations, we get 

f(b{T,p)<^-^) ^ expi-t) for T = 256 ( --^ + ^ + ^-^) + ( + S/]) t 

□ 

The next result still concerns the "dynamic" expansion of any subset of nodes. It will be 
applied when the subset of informed node is large (say at least n/2) in order to get a good 
bound on the completion time of the last phase of the flooding process. Let us define the 
following r.v. 

1 if 3t' 3j e ^ : r < r' ^ r + r and {i, j} G E^' , . 

otherwise 

Lemma 12 If Q{[n]^{Et}t^Q) is {M, a, (3) -stationary then, Vr, t 1 it holds that, for every 
A C [n] and for every i € [n] \ A, it holds that 

P(eS = o|^^r) ^ exp(-t), where T = 2 + t 

Proof. For brevity's sake, we omit the conditioning by E^r- For every s = 1, . . . ,T, define r.v. 

Observe that Yg = fs{E^r+s) for a suitable function fg- Lemma [9] implies 

\A\a 



P (n = 1 I Er+s-i) = P (deg[+^ > I ^ p 



2 + 2\A\ap 

By applying Lemma [7] to E-^+i, . . . , Et-+t and Yi, . . . ,Yt, we get 

P (^^i; = 0^ ^ P {B{T,p) = 0) = {l-pf ^ exp(pr) = exp(-t) 



□ 



C.l Flooding in Dynamic Graphs 

In what follows, we bound the time required to obtain at least n/2 informed nodes. This is 
the spreading phase. 
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Lemma 13 (Spreading Phase) IfQ{[n],{Et}t^o) is {M , a, (3) -stationary then, \/t ^ T with 
T = O (^(^ + fSflog^ it holds that 

^(I'-l < 2) < S5 

Proof. (Sketch of). Observe that, for any |^|, the bound on T for t = B(logn) in Lemma [TT] 
satisfies 

\\|>i|n^a^ na n J \\A\na J J \\na J I 

From Lemma [TH after every time interval of T epochs, the size of the set of informed nodes at 
least doubles as far it is smaller than n/2, with high probability (say 1 — 1/n'^). By a simple 
application of the Union Bound, we get that, after a O(logn) number of such time intervals, 
with high probability the number of informed nodes is at least n/2. □ 

Lemma 14 (Saturation Phase) Let Q{[n], {Et}^^^) be {M , a, (3) -stationary and assume the 
flooding process is in some epoch r such that ^ n/2. Then, w.h.p. all nodes get informed 
within O ((^ + /3) logn) epochs. 

Proof. (Sketch of). Assume that we are in some epoch r where |It-| ^ n/2, then by choosing 
A = Ir and t = 0(logn), Lemma [T2] implies that, with high probability, every node gets in- 
formed within O ((^ + /?) logn) epochs. Then, by the Union Bound, all nodes gets informed 
after the same number of epochs. 

n 

Proof of Theorem [H Thanks to Lemma [T3l w.h.p., after O (^M (^ + log^ n^ steps the 

set of informed nodes will be at least n/2. Then, from Lemma [TH we have that, after further 
O (M + /?) logn) steps, w.h.p. all nodes will get informed. 

□ 



D Flooding in Node-MEG: Proof of Theorem [3] 

We firstly need some notations. For every x ^ S, define 

T{x) = {y(^S\C{x,y) = l} 

In words, V{x) is the set of states that are at one hop from state x. For any node z £ [n], let vi 
be a probability distribution over the set of states S. The symbol v (without subscript) will 
denote the product probability distribution over S. Assuming that nodes are random 

with the probability distribution u, for every i,j € n] and for every A G [n] — {i}, define 
binary random variables e^j and e^^, 

e^ j = 1 if nodes i and j are connected 

ej'^ = 1 if node i is connected to some node in A 
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It is easy to verify that the followings hold 



x€3 

where, for any binary random variable X, X is the complementary random variable defined 
as X = 1 iff X = 0. Notice that VjiT{x)) is the probability that node j is connected to a 
fixed node in state x. Let vr be the stationary probability distribution of the Markov chain 
M.. With an abuse of notation, we denote by tt also the probability distribution over nie[n] ^ 
given by the product of vrs. For the sake of simplicity, we omit the explicit dependence on tt 
of random variables and probabilities. Thus, we write Cjj, Cj^^ to mean, respectively ejp ejj^. 

Lemma 15 Let NM = NM{n,M,C) be a node-MEG such that 

-Pnm2 ^ r] (Pnm)^ for some i] ^ I 

Then, 

Vi, j G [n]V^ C [n] - {i, j} P {e^^A ■ e,,^) ^ 17r/P (e^.^) P (e^, a) 

Proof. For the sake of convenience, let q{x) and q{x,y) denote, respectively, Tr{T{x)) and 
7r(r(x) ur(y)). Define 

V= \ xeS\ q{x) > ^ 



1^1 J 

Claim 1 For every k & [n] — A, it holds that 

P (eM) ^^(pnm-Y1 <^)Qi^)] + 
\ xev / 

Proof. We upper bound P {ej^). Clearly, it holds that 

P(MA) = Y.7r{x){l-qix))\^\+ ^ 7r(x)(l - (8) 
xev xes-v 

li x eV, then 

(1 - ^ e-«(^)l^l ^ e"^ ^ e"^ (9) 

lixeS-V, then 



(1 - qix)r\ ^ e-«(^)l^l < 1 - ^ = 1 - ^q{x) (10) 

2^y\A\ ^ 
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where we used the inequahty e ^ ^ 1 — x/{2a) (that holds for q ^ 1 and 
^ X ^ a). By combining InequaHties dH), ([9]), and (fTO|) . we obtain 

P (eiM) < ^ E ^(^) + E vr(x) ( 1 - V^g(x) J 

= ^^(^) + E ^(^) ( 1 - ^«(^) ) - E ^(^) ( 1 - 



1^1 , , 



ivr(F) + 1 - ^Pnm - E ^(^) + ^ E 



^ 1 - l^V) - ^ ( Pnm - E 



and the thesis immediately follows. 
We distinguish two cases. First we assume that the following holds 



2' NM 



It holds that 



— — — - ( 7r{x)q{x) j ^ Tr{x)q{x)'^ (by Jensen's inequality) 
^ \xev J xev 

^ V (-Pnm)^ (by lemma's hypothesis) 

< 4:7] TT{x)q{x) (by Assumption ([TT]) ) 

\xev J 



From this, it is immediate to derive that 
Thus, thanks to Claim [1] we obtain 



P (efc,A) ^ ^ (^Pnm - E + ^^(^) ^ ^ 



It follows that 

P {ei^A ■ ej,A) ^ P {e^,A) ^ ^ 877? (e,,^) P (ej, a) 
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Now, we consider the opposite case: 

J2 7rix)q{x) ^ ^Pnm (12) 



xev 



For any two binary r. v. X and Y, it holds that X Y = {X \/ Y) and thus P{X -Y) = 
P (X • F) + 1 - P (X) - P (F) . Hence, 

P {ei,A ■ ej,A) = P {eij ■ ejj) + 1 - P (e^) - P (e^) (13) 
First, we focus on upper bounding P (e^^ • e^). It holds that 

V{e^-ej;A) = Y.Y. <^)<y) n ^(-5 - (r(a^) U T{y))) 
xeSyeS heA 

= J2 7r(ar)7r(y)(l-7r(r(x)ur(y)))l^l 

x,yeS 

Define 

R = {{x,y)eSxS\T{x)nT{y)^il}} and R = SxS-R 
Observe that if {x,y) G R then q{x,y) = q{x) + q{y)- Thus, it holds that 

P{eU-eJj)= vr(x)7r(y)(l-(z(x,2/))l^l+ Yl vr(x)7r(y)(l - g(x, y))!^! 

= Y <xMy){l-q{x)-q{y))\''\+ Yl Ax)n{y){l - q{x,y))\^\ 
{x,y)eR {x,y)eR 

= Ai + A2 (14) 

where 

ki=Y 7r(a;)7r(y)(max{l - q{x) - q{y), 0})l^l 

x,y(^S 

A2= Y vr(x)7r(y)[(l-g(x,2/))l^l-(max{l-g(a;)-g(2/),0})l^l 

{x,y)eR 

To bound Ai we use the inequality max{l — a — 6,0} ^ (1 — a) (1 — 6), that holds for any a 
and b with ^ a, 6 ^ 1: 

Ai^ ^ 7r(x)7r(y)(l-g(x))l^l(l-g(y))l^l 

x,y&S 

a;e5 yeS 

= P(ei:A)P(^) (15) 
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By combining Inequalities (fT3|) . p^ . and ([IS]), we obtain 



P iei,A ■ ej^A) ^ 1 - P (ei;i) - P (e^) + P (elj) P (e^^) + A2 
= (1 - P (elj)) (1 - P (e77)) + A2 
= P {e^,A) P {ej,A) + A2 

To upper bound A2, we use the following 

Claim 2 // ^ ^ 1 and 1^1 — a + then, for any integer k ^ 1, 

(l-a + b)'' - (max{l - a, O})'^ ^ kb 

Proof. By distinguishing the two cases 1 — a ^ and 1 — a > 0, and in the latter 
by induction on k. □ 

Let q{x, y) denote 7r(r(x) n r(y)). Observe that 

q{x, y) = q{x) + q{y) - q{x, y) 
From Claim [2] with a = q{x) + q{y) and h = q{x,y), we get 

(1 - q{x, - (max{l - q{x) - q{y), 0})l^l ^ \A\q{x, y) 

It follows that 

A2 ^ |A| ^ 'K{x)'K{y)q{x,y) 

A\ ^ 7r(x)7r(y)g(x, y) (since (x, y) ^ R ^ q{x, y) = 0)) 

x,yeS 

A\ ^ 7r(x)7r(y)^7r(z)C(z,x)C(z,y) 

265 2:65' yeS 

A\Y,<z)q{zf 

zeS 

^ rj\A\ (Pnm)^ (from lemma's hypothesis) 



By combining this with Inequality (jl6p . we obtain 



P {e^,A ■ ej,A) ^ P {e^,A) P (ej-,A) + (^i 



From Claim [T] and Hypothesis (|12p we get 

P (efc.A) ^ ( i=^M - E ^(^)^(^) ) ^ 



NM) 
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In conclusion, from this and Inequality (jl7p we get the thesis: 



P (ei,A • ej, a) ^ P {ei,A) P (cj-.a) + r]\A\ {Pi 



NMJ 



^ P (e,,A) P (ej,A) + il\A\—j= 

= P (ei,A) P (ej,A) + 16r?P (ei,A) P (c^.a) ^ 17r/P (e,,A) P (e^, a) 



□ 



Lemma 16 For every i = 1, . . . ,k, let tpi and Q be any two probability distributions over any 
domain Let tp and ( denote the product probability distributions over Y[i=i of tpiS and 
Qs, respectively. Then, it holds that 

k 

\\i^-C\Wy^^U^-C^\\TY 
i=l 

Proof. Let Q, = Y[i=i ^i- We denote (xi, . . . , Xk) by x. For every i, let = 11^=1 j^i ^^^'^ 
let x~* denote (xi, . . . Xj+i, . . . ,Xk). It holds that 

IIV'-CIItv = ^EI^(^)-C(^)I 

=^E nv'.(^.)-no(x,) 

xen j=i 3=1 

^ fc k i—1 ,„ ,„ . J. 

= 2 E n ^^-(^^o + E n n - n o(^i) -eh n ^^■(^^■) 

xef2 j=l i=2 j=l j=i j=l 

=^E Eno(^.)nv',(^.)-Erio(^.) n ^^^^^ 



k i-l 



i=2 j=l 



xGQ i=l j=l j=i 
k i-l 



i=l j=l j=i+l 
k k i—1 



^ E E ^» ) n ) n (a^j ) - X] Ci ) n o ) n ^/'^ ) 

x£fl i=l j = l j=i+l i=l j = l j=i+l 

k i—lk 



xgf2 i=l j=l j=i+l 

^ fc i—1 k 

^ 2 E E ~ n n ^j(^j) 

x£U i=l j=l j=i+l 

=^E E \Mx^)-ux.)\ no(^.) n ^^■( 



(18) 
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Observe that 

i— 1 k i—1 k 

E no(^.) n ^^■(^^•)= E no(^.) e wu^^) 

E no(^.) 

(xi,...x,_i)en-=lf^j^=i 

= 1 

Hence, from this and InequaUty psp the thesis fohows. □ 
Lemma 17 Zei NM = NM{n,M,C) he a node-MEG such that 

-Pnm2 ^ V [P^mf for some r/ ^ 1 (19) 
Moreover, let v he a product probahility distribution such that 

w ^ r 1 II II ^ (-Pnm)^ 

V? e n - ^ TV ^ — (20) 

2n 

Then, for any i,j £ [n] and for any A C [n] — {i,j}, 
(^) P ^ ^ 

Proof. Let Cjj = {(^i, . . . , x„) G HILi ^ I C'(^i) ^i) = !}• It holds that 
P (er,,) = i^(Q,,) 

^ 7r(Ci ,) - n ^^^^^ (from ([20]) and LemmalHl) 
2n 

n (-Pnm)^ ^ -Pnm 
= Pnm ^ ^ ^ 

This proves (jaj). 

For any /c G [n] and for any A C [n] — {A;}, let Ck,A = {(a^ij ■ ■ ■ ) Xn) € nr=i 5* | 3/i G ^ : 
C{xk,Xh) = 1}. For any i,j G [n] and for any A O [n] - {i,j}, let Cij,^ = {(xi,... ,x„) G 
nr=i S \3h,k e A: C{xi,Xh) = 1 A C(xj, x^) = 1}. It holds that 

^ vr(Ci a) + n ^'^™^ (from mh and Lemma [ISD 
2n 



(-Pnm 



i2 



P(ei,^-e,-A) + ^^r^ (21) 
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Since Hypothesis (jl9p holds, Lemma [T5] ensures that 

P {ei^A ■ ej^A) ^ 17r/P (gj,^) P (ej,^) 
From this and Inequahty (j2ip we obtain 

P «A • ^ 17r?P (e,A) P (e,- a) + 

^ 18r/P (ei,A) P (ej,A) (since Pnm ^ P (ei,A) , P iej,A)) 

= 18r,7T{Ci,A)7T{Cj,A) 

^ 18r,{u{a,A) + ^^^^){HCj,A) + ^^^2^ ) (from (I2D|) and Lemma^ 

^ {elA) + (eU) + ^) (22) 

Since Pnm = ^ i^ij) = 7r({(2;i, . . . , x„) G nr=i I C{xi,Xj) = 1}), from (j20]) and Lemma [TBI 
it derives that Pnm ^ P (ej'j) + ((Pnm)^)/2. It follows that Pnm/2 < Pnm - ((Pnm)^)/2 < 
P (^e'^j^ ^ P {^iAj ■ Fi'oni this and Inequality ([2^ we get 

P «A • ^ 18r?(2P (er,A))(2P (ej,^)) = 72r?P (ej;^) P (ej,^) 

□ 

Proof of Theorem [S], 

Let Q = G{[n],{Et}t^o) be the dynamic graph yielded by NM{n,Ai,C) and let e*^- and e*^ 
be the relative random variables. Consider M = Tmix log(2n/(PNM) ) be the duration of an 
epoch. Prom standard results on Markov chains, it derives that, for any node i, whatever the 
probability distribution at time t the probability distribution i/i at time t + M is such that 

Ui,,- ^ 2-i°s(2n/(PNM)^) = (-Pnm)' .^S) 

2n 

This and the theorem's hypothesis satisfy the hypotheses of Lemma [T71 and in turn it implies 
that 

P(el,f I i?,^(._,)M)^^ (24) 

P {e^J ■ e]^i I Pt^(._i)M) ^ 7277P {el^ \ Pi^(._i)M) P [efi \ Et^ir-m) (25) 

Hence, Inequalities ([25|) and (plj) show that ^ is a (M, Pnm/2, 7277)-stationary dynamic graph. 
Then, from Theorem [H the flooding time in Q is, with high probability, 

2\ ^ , vo^ \ i„„2 



log(2n/(i^M)') + 72r/ log' n 

and, taking into account the hypothesis Pnm ^ l/n*-*(^\ it is 

O (Tmix f^— + log^n 
\ V^Pnm / ; 



□ 
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D.l Flooding in Classic Mobility Models: Proofs 
Proof of Corollary [4l 

Let vr be the stationary distribution of the Markov chain A4 of NM. For any x € S, let u{x) 
be the point in TZ where a node is when its state is x. By assuming that the states of the 
nodes are random with vr, let q{x) be the probability that a fixed node is connected to another 
fixed node being in state x. Clearly, 

yes : u(y)eD{u{x),r) 

Since NM is a sufficiently refined discrete version of the random trip model T, it holds that 



q{x) ^ / Fr{u)du (26) 

J D{u(x),r) 

From Hypothesis (jb]), for every v G Br, 
This implies that 

vol(i:'(7;,r)) Cdr'^ 

where q is a constant depending only on d. By combining (j26|) and (j27p . we have that, for 
every x £ S with u{x) G Br, 



It holds that 



-Pnm = y^7r(x)g(x) 

^ ^ 7r(2;)g(x) 

> ^^^lijl^ Y ^(^) Inequality ([M])) 



xes : u(x)eBr 

d 

F-Y{u)du 



(5vol(7^) Jrs 

Cdr'^ vol (Br) u +1. • /iKh 

^ (from Hypothesis (|bD 



(5vol(7^) (5vol(7^) 



> X2^^77m Hypothesis ©) (29) 

o^voi(/<.j 
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From (j26p and Hypothesis ((aj), it derives that, for every x £ S, 

qix) ^ I Fr{u)du < = ^ (30) 



D{u{x),r) 



It holds that 



-Pnm2 = y^7r(x)g(, 

x&S 



xf 



( 6C(ir \ 

~ ( — I / tt(x) (from Inequahty (130])) 



~ vol(7^)2 

From this and Inequahty (j29|) . we get 



^^^^ - ^;d(^ - V^^^^y -A^^^^^^ 

It follows that the node-MEG NM satisfies the hypotheses of Theorem [3] with r/ = 6^ /}? and 
thus, with high probability, the flooding time is 



where we have used Inequality (129p . □ 
Proof of Corollary O 

It is easy to see that, any node-MEG yielded by a random-path model is a Markov Trace 
Model (MTM), a general class of models introduced in [T3]. Since RP is simple and reversible, 
Theorem 11 in [14j implies that the stationary distribution vr of is uniform. For any state 
x G S*, let u{x) G y be the point where a node is when its state is x. By assuming that 
the states of the nodes are random with vr, let q{x) be the probability that a fixed node is 
connected to another fixed node being in state x. Since vr is uniform and RP is simple, it 
holds that 

q{x)= <^) = ^\\{y^S\u{y)=u{x)}\ = tE^^ (31) 

It follows that 

Pnm = = ^ E ^-pHx)) = tLY. #^(^)' (32) 

where the last equality derives from \{y £ S \ u{y) = n}| = ^-p^u). Thanks to Jensen's 
inequality it holds that 

Euev#v{uf ^ fEueV#Pi^)V _ \S\' 



v\ V 11^1 / w 
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Thus, from ([32]) . we have 

Pnm>t^^ (33) 



From dST]), it holds that 

pnm2 = Yl ^i^u^f = ^ E #^("(^))' = ^ E #^(^)' 

Since RP is (5-regular and ()33p holds, it follows that 



From this and taking into account the hypothesis \V\ ^ n^^^^ and (j33p . Theorem [3] can be 
applied with rj = 6'^ obtaining that, with high probability, the flooding time is 



O ( T^ix( ^+<5^ 1 log^n 



□ 
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